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This  work  contains  the  results  of  a study  aimed  at  the 
development  of  two-  and  three-dimensional  numerical  proce- 
dures for  computing  the  flowfield  generated  by  the  inter- 
action of  a blast  wave  and  a rigid  body.  A number  of  numer- 
ical procedures  were  applied  to  two-dimensional  problems 
including  both  implicit  and  explicit  algorithms.  Each  was 
tried  on  the  blast  wave-cylinder  interaction  problem. 

MacCormack’s  method  with  added  fourth-order  dissipation 
yielded  the  best  results  and  was  then  applied  to  the  blast 
wave-truck  interaction  problems  in  two  dimensions. 

MacCorraack's  method  was  also  used  in  three  dimensions  to 
determine  the  flowfield  that  results  when  a blast  wave 
strikes  a rectangular  parallelepiped  at  an  arbitrary 
angle.  Both  the  two-  and  three-dimensional  computations 
were  compared  with  experiments  in  a number  of  ways.  Two 
dimensional  density  contours  show  qualitative  agreement 
for  shock  front  location  and  Mach  stem  fomation  with 
spark  shadowgraphs  taken  in  a shock  tube. ^Pressure-time 
histories  indicate  good  quantitative  agreement  between 
theory  and  experiment  both  in  two-  and  three-dimensions. 

INTRODUCTION 

The  accurate  prediction  of  the  effects  of  blast  waves  encountering  vehicles 
and  structures  is  essential  in  the  design,  survivability,  and  hence  effective- 
ness of  these  configurations.  Detailed  experimental  blast  wave  interaction 
data  is  both  costly  and  difficult  to  obtain.  Moreover,  these  experiments 
frequently  do  not  provide  a complete  picture  of  the  blast  wave  interaction  flow- 
field.  Actual  experiments,  in  fact,  only  yield  pressure  data  at  a few  selected 
points  on  the  models.  As  a consequence  essential  design  parameters  are  often 
difficult  to  define. 

Design  information  for  the  kind  of  blast  wave-vehicle  encounter  as  pictured 
In  Fig  1 may  be  obtained  experimentally  in  several  ways.  In  one  approach  a 
targe  explosive  charge  is  detonated  near  on  instrumented  vehicle  or  structure 
enabling  direct  measurements  to  bo  taken.  This  approach,  while  realistic,  i* 
costly,  results  in  few  data  points,  and  frequently  provides  little  understanding 
of  the  important  flowfield  phenomena.  Another  approach  is  to  place  a model  of 
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the  target  inside  a shock  tube.  In  this  case,  better  control  can  be  provided 
and  the  experimental  cost  is  less.  Nevertheless,  the  experiments  are  still 
limited  in  their  range  of  applicability  as  a result  of  scaling,  shock  tube  wall 
effects,  realistic  wave  shapes  and  flow  duration. 

As  an  alternative  to  the  experimental  description  of  the  blast  wave  inter- 
action phenomenon,  one  can  use  computational  fluid  dynamics.  This  is  the 
approach  adopted  here.  Accurate  finite  difference  simulations  offer  the  possi- 
bility of  providing  design  data  at  a relatively  low  cost.  Such  a simulation 
provides  a complete  flowfield  description  that  is  essential  to  a fundamental 
understanding  of  the  fluid  mechanics  and  a necessity  for  an  effective  structural 
design.  The  numerically  generated  flowfield  data  can  then  be  integrated  to 
yield  other  vital  information  such  as  the  total  loads,  center  of  pressure,  and 
overturning  moments. 

In  the  past  second-order  finite-difference  procedures  have  been  used  by 
Kutler,  et  al.  (See  References  1 and  2)  to  solve,  simple  shock-diffraction 
problems  involving  both  regular  and  Mach  reflections  of  the  Incident  shock.  In 
these  cases  all  discontinuities  wore  fit,  i.e.,  treated  as  sharp  discontinuities. 
Blast  wave  encounter  problems  with  two-dimensional  wedges  and  three-dimensional 
cones  in  supersonic  flight  have  boon  solved  by  Kutlor,  ot  al.  (See  References  3 
and  4)  again  using  second-order  finite -difference  procedures.  In  those  instances 
a "shock-capturing"  philosophy  was  employed  and  resulted  in  an  accurate  descrip- 
tion of  the  so-cailed  "shock-on-shock"  problem. 

lr.  the  present  paper,  these  "shoek-eaptuvung"  flowfield  simulation  tech- 
niques have  been  adapted  to  the  blast-wave  interaction  problem.  Only  inviscid 
flow  problems  have  been  considered,  but  both  complex  two-dimensional  oral  simple 
three-dimensional  geometric  configurations  have  been  used  as  targets, 

GOVERNING  EQUATIONS 

Several  assumptions  are  made  in  the  present  study  of  hirst  wave  encounters 
with  targets.  The  first  is  that  the  blast  wave  Is  assumed  to  he  planar  relative 
to  the  target  and  that  conditions  behind  the  wave  can  bo  adequately  and  consis- 
tently described.  Secondly,  viscous  effects  are  ignored,  finally,  any  effects 
which  result  from  radiative  heating  on  the  target  are  assumed  negligible,  and  a 
perfect  gar  equation  of  state  is  employed. 

Under  the  above  assumptions,  the  governing  partial  differential  equations 
are  the  unsteady  Euler  equations.  To  permit  the  mapping  of  two-  or  three- 
dimensional  complicated  physical  regions  into  rectangular  or  cubical  computa- 
tional douains  respectively,  the  following  independent  variable  transformation 
is  employed: 

t « t;  C - £(t.  x,  y.  i);  n - n(t,  x,  y,  s);  c * t(t,  x,  y,  a)  (1) 

Because  the  above  transformation  maps  the  body  and  outer  boundary  surfaces  onto 
constant  coordinate  lines  and  planes,  application  of  the  boundary  condition 
procedures  is  facilitated.  The  above  transformation  also  permits  (he  clustering 
of  grid  points  in  the  vicinity  of  the  body. 

Under  this  transformation,  the  governing  partial  differential  equations  in 
strong  conservation  law  fora  become 

\ 4 % * Gt  * 0 

where  the  flux  vectors  t,  f,  G of  liq.  (2)  assume  the  form 


(2) 
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> w£z  , V = nt  + ur^  4'VHy  wnz  , and 

are  the  contravariant  velocity  components  without 


In  the  conserved  variables  of  Eq.  (2)  p represents  the  pressure,  p the 
density,  u,  v,  and  w the  three  Cartosian  volocity  components,  and  e the  total 
energy  per  unit  volume.  For  an  ideal  gas,  the  pressure,  density,  and  velocity 
components  are  related  to  the  energy  by  the  following  equation: 


e 


+ o 


v2* 


(4) 


The  metrics  required  by  Kq.  (3)  in  general  are  not  known  analytically  and 
must  be  evaluated  numerically.  Details  for  this  procedure  are  given  in 
Reference  S, 


NUMERICAL  ALGORITHMS 

The  transformed  governing  equations  (Eq.  (2))  were  solved  by  both  explicit 
and  implicit  f inite-difference  procedures.  These  schemes  included  NacCormaek's* 
explicit  method  with  an  additional  fourth-order  dissipation  term,  Beam  and 
Warming’s f> 8 implicit  method  in  the  delta-form,  and  Steger  and  Warming's’ 
explicit  upwind  scheme.  MaeCoraack's  scheme  captured  the  shock  within  the 
least  number  of  grid  points  and  consumed  the  least  amount  of  machine  time,  and 
therefore,  it  is  the  only  procedure  presented  here,  the  others  can  be  found  in 
Reference  S. 

MacCoraack's  method  is  a second-order,  noncentcred  predictor-corrector 
scheme  and  appears  as  follows: 


q * q”  - At  (Aj,En  ♦ A^p"  * &^) 

($) 

qn*1  - |{q  • q”  - AtfC^E  * V^F  ♦ V^S)  ♦ e.Dn] 

where  ft  implies  that  the  flux  vector  E is  evaluated  using  elements  of  the 
predicted  value  q,  and  A and  v are  the  standard  forward  and  backward  difference 
operators,  the  quantity  I)  represents  a fourth-order  dissipation  term  in  all 
thre-  directions  whose  effect  is  governed  by  the  dissipation  constant  c. 

GRID  GENERATION 

The  generalised  coordinate  transformation  given  by  Eq.  (1)  permits  the  use 
of  grids  based  not  only  on  standard  coordinate  systems  such  as  cylindrical  or 
spherical  hut  also  numerically  generated  grids  such  as  those  obtained  by  solving 
elliptic  partial  differential  equations.  In  this  study  both  analytically  and 
numerically  detenained  grids  were  used  to  discretiac  the  physical  regions  of 
interest . 
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For  two-dimensional  targets  consisting  of  cylinders  and  rectangles,  an 
analytically  described  mesh  based  on  a cylindrical-like  coordinate  system  is 
employed  (see  Figs.  2a  and  2b).  Points  can  be  clustered  near  the  body  for  better 
resolution  by  using  an  exponential  function. 

In  order  to  discretize  an  arbitrary  two-dimensional  shape  such  as  a truck 
and  its  surrounding  flowfield,  a numerically  generated  grid  was  used  (see 
Fig.  2c).  The  partial  differential  equations  used  for  this  process  are  elliptic 
and  satisfy  a maximum  principle.  The  technique  for  generating  these  grids  was 
obtained  fTom  R.  L.  Sorenson  at  NASA's  Ames  Research  Center  and  was  a result  of 
a paper  by  Steger  and  Sorenson.10  Details  of  the  mesh  generation  procedure  can 
be  found  in  that  paper,  but  suffice  it  to  say  that  the  method  can  treat  arbitrary 
bodies  and  outer  boundaries  with  clustering  near  the  body. 

For  the  three-dimensional  problem,  the  flowfield  resulting  from  the  inter- 
action of  a blast  wave  and  a rectangular  parallelepiped  was  desirod.  Discreti- 
zation of  this  flowfield  was  accomplished  using  a spherical -typo  grid  as  shown 
in  Fig.  3. 

BOUNDARY  AND  INITIAL  CONDITIONS 

Two  types  of  computational  boundaries  were  considered;  those  across  which 
there  is  no  mass  flow  (impermeable  surfaces)  and  those  across  which  there  is 
mass  flow  (permeable  surfaces).  Impermeable  boundaries  include  solid  walls, 
planes  of  symetry,  and  slip  surfaces  whereas  permeable  boundaries  include  shock 
waves,  porous  walls,  and  inflow  or  outflow  boundaries. 

In  the  blast  wave  interaction  problems  of  interest  here,  boundary  condition 
procedures  are  required  at  the  body,  the  planes  of  symmetry,  and  along  the  outer 
boundary.  For  inviseid  flows,  the  boundary  condition  at  the  surface  of  a body 
requires  flow  to  be  tangent  to  the  body.  This  implies  that  the  velocity  compo- 
nent V of  Rq . (S)  must  be  equal  to  zero  at  the  body.  In  order  to  simulate  this 
numerically,  the  image  plane  concept  is  used.  By  employing  this  concept,  an 
image  line  of  nodal  points  is  established  which  falls  one  mesh  Interval  inside 
of  the  body  (see  Pigs.  213).  Flow  variables  along  this  line  are  obtained  by 
use  of  the  flow  variables  at  the  body  and  the  adjacent  flowfield  interior  points, 
both  at  the  previous  time  step.  The  flow  variables  for  the  new  time  step  at  the 
body  ean  now  be  obtained  by  the  same  numerical  algorithm  that  was  used  for  the 
interior  points,  ami  hence  with  the  use  of  the  implicit  procedure,  the  body 
points  are  also  updated  implicitly.  The  details  of  this  procedure  for  both  the 
two-  ami  three-dimensional  cases  are  given  in  Reference  S, 

To  numerical ly  simulate  the  ground  plane  (plane  of  symmetry)  in  both  the 
two-  and  three-dimensional  problems,  the  reflection  principle  is  employed.  In 
this  approach  the  pressure,  density,  energy,  and  tangential  velocity  components, 
are  treated  as  even  functions  with  respect  to  the  ground  plane  while  the  normal 
velocity  component  is  treated  as  an  odd  function, 

RliSUt.TS  AND  DISCUSSIONS 

In  this  section  numerical  results  are  presented  which  describe  the  inter- 
action of  a blast  wave  with  a body.  In  tvo-dimortsions,  a comprehensive  set  of 
results  is  presented  for  the  blast  wave-cylinder  interaction  while  partial 
results  arc  presented  for  the  blast  wave-truck  interaction.  In  three-dimensions 
partial  results  are  also  presented  for  the  bias*  wave  parallelepiped  inter- 
action problem. 

Numerical  and  experimental  results  for  the  blast  wave-cylinder  interaction 
problem  arc  presented  in  severs 1 ways.  In  Fig.  *a  computed  pressure-time  his- 
tories are  compared  with  experimental  result*  obtained  by  Pearson  ct  al,s  at  the 
Ballistic  Research  laboratory.  The  inset  in  each  figure  depicts  the  measuring 
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station  of  the  transducer.  The  blast  wave  passes  over  the  cylinder  from  left  to 
right.  In  each  case,  the  smooth  curve  is  the  computed  result  while  the  oscilla- 
tory curve  is  the  measurement.  The  oscillations  represent  excursions  from  true 
values  due  to  transducer  ringing,  local  unsteadiness  due  to  model  oscillations, 
general  flow  anomalies,  etc.  A sample  of  the  more  interesting  cases  with  signif- 
icant wave  structure  are  shown  (8  = 45°,  135°,  180°).  All  the  locations  (except 
for  180°)  have  a similar  structure.  A sharp  rise  in  the  pressure  signifies  the 
arrival  of  the  incident  blast  wave  at  the  measuring  station.  This  is  followed  by 
an  exponential -like  decay  and  a second,  more  "smeared"  jump  to  a plateau  value. 
The  delays  from  t = 0 on  the  abscissa  represent  the  time  before  the  shock  arrived 
at  the  sampling  station.  The  results  can  be  more  fully  explained  with  the  help 
of  Figs.  4b,  computed  isopycnics,  and  4c,  spark  shadowgraphs  (obtained  by 
Vandromme'2) . Interesting  results  develop  as  the  shock  encounters  the  cylinder 
and  passes  over  it.  A triple  point  and  Mach  Stem  appoar  partway  up  the  windward 
side  of  the  cylinder  as  is  clearly  seen  in  both  the  computed  isopycnics  and  the 
shadowgraphs.  The  contact  surface  emanating  from  the  triple  point  (as  seen  in 
Fig.  4c)  is  not  reproduced  in  the  isopyenic  plots  of  4b.  This  is  because  of  in- 
adequate resolution  of  the  computational  grid.  Better  results  would  be  obtained 
with  an  increase  in  the  number  of  grid  points  or  an  adaptive  gridding  scheme. 

As  tho  primary  shock  continues  toward  the  rear  of  the  cylinder,  it  eventually 
reflects  and  is  propagated  upstream.  As  this  reflected  shock  arrives  at  a mea- 
suring station,  it  gives  rise  to  the  second  "smeared"  jump  in  Fig,  4a.  The 
smeared  nature  of  tho  jump  is  manifest  in  viscous  effect  near  the  surface  which 
now  begin  to  make  their  effects  felt.  The  discrepancy  in  the  pressure  amplitude 
in  the  bottom  figure  of  4a  and  the  latter  half  of  the  top  two  figures  is  also 
attributable  to  viscous  effects.  These  effects  were  neglected  in  the  computa- 
t ion . 


The  interaction  of  a blast  wave  with  a truck  for  intermediate  times  is 
shown  in  Fig.  5 in  the  form  of  pressure  contours.  For  this  ease  the  grid  site 
consisted  of  52  points  in  the  j-direction  (along  the  body)  and  25  points  in  the 
k-direct.ion  (normal  to  the  body).  The  blast  wave  had  a strength  of  34. S kPa 
overpressure.  Figure  Sa  shows  the  reflected  blast  wave  free  the  front  of  the 
truck  while  Fig.  Sb  shows  a reflected  wave  from  the  windshield.  In  Fig.  Sc  the 
blast  wave  is  beginning  to  expand  over  the  cab  of  the  truck  while  the  reflected 
waves  from  the  front  and  windshield  move  away  from  the  vehicle.  This  particular 
calculation  does  demonstrate  the  versatility  of  the  arbitrary  body  mesh  gener^ 
ator  amt  shock -capturing  ability  of  MacConiack's  method  for  computing  compli- 
cated flowfields. 

The  three-dimensional  interaction  of  a blast-wave  with  a rectangular  paral- 
lelepiped was  computed  using  KacCoroack's  method  with  fourth-order  smoothing. 

For  this  calculation  the  grid  consisted  of  S4  points  in  the  j-direction  or 
around  the  body,  24  points  in  the  k-direction,  and  16  (mints  in  the  {-direction, 
or  between  the  body  and  the  outer  boundary  (see  Fig,  5). 

Numerical  results  were  obtained  for  which  experimental  data  was  available. 
The  experimental  data  was  obtained  by  placing  pressure  gages  at  various  posi- 
tions on  the  model  and  recording  the  time  histories  of  the  pressure  as  it  was 
struck  by  the  blast  wave  (see  Fig.  5). 

The  Initial  conditions  for  this  case  consisted  of  a blast  wave  Mach  number 
of  1.14,  freostreaa  pressure  of  101.33  kPa,  and  angle  of  incidence  a of  52.5*. 
the  distance  r^  from  the  origin  to  the  blast  .wave  fscc  Fig.  3)  was  4.91. 
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The  results  of  this  first  calculation  are  shown  in  Fig.  6 where  the  pressure 
in  kilo-Pascal s is  plotted  as  a function  of  time  in  milliseconds  for  three  sta- 
tions. The  solid  dots  shown  on  each  of  the  curves  in  Fig.  6 is  the  experimental 
data.  The  agreement  is  acceptable  and  can  be  made  much  better  with  better  grid 
resolution. 
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